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Abstract 

We calculate the Lamb, fine and hyperfine shifts in bb with n = 2, 1 = 0, 1. 
Radiative corrections as well as leading nonperturbative corrections (known 
to be due to the gluon condensate) are taken into account. The calculation is 
parameter- free, as we take A, (a s G 2 ) from independent sources. Agreement 
with experiment is found at the expected level ~ 30%. Particularly inter- 
esting is a prediction for the hyperfine splitting, M average (2 3 P) — M(2 1 P 1 ) = 
1.7 ± 0.9 MeV, opposite in sign to the cc one (« — 0.9MeV), and where the 
nonzero value of (a s G 2 ) plays a leading role. 
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1 Introduction 



In a previous paper |l] (hereafter to be referred as TYq) we presented an 
evaluation of the potential for heavy qq systems fl], |2|. The evaluation in- 
cluded relativistic effects, one-loop radiative corrections and (for the spin- 
independent part) the dominating two-loop ones. With this we evaluated 
a number of quantities, taking into account also leading nonperturbative 
corrections, which are known to be due to the contributions of the gluon 
condensate. It was shown that a very good account could be given of the low- 
est lying bb bound states (some features of cc were also discussed). Notably, 
both the energy and wave function (this last through e + e~ decay) of the 
states with n = 1 were given; the splittings between these states and those 
with n = 2, / = 0, 1 were reproduced in what is essentially a zero parameter 
calculation using only the known values of the basic QCD parameters, 

A{n f = 4, 2 loops) = 200 Ijjg MeV 

(a s G 2 ) = 0.042 ± 0.020 GeV 4 (1) 
m b = 4906 ±j» ~\ t\l MeV. 

Actually we preferred in TY to deduce m\> from the mass of the T(15') state. 
The errors given for this quantity in ([]]) correspond to that in A (the first), 
to that in the gluon condensate (the second); the third is an estimated sys- 
tematic error. 

The value of given in (Jl|) is for the 'pole mass, which is the appropriate 

quantity to be used in a Schrodinger equation. It corresponds to a running 
*We will freely use the notation of TY 
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mass value of 

m b (m 2 b )= 4397 + 7 2 ~l t^MeV, (2) 

which compares favorably with the SVZ estimate jl) of 4250 ± 100 MeV. 

For some of the states with n — 2, I — 1, no result could be given; only 
the perturbative contributions were presented and they failed to reproduce 
the experimental values. This was because the nonperturbative corrections, 
more involved than for the n = 1 case, had not been calculated at the time. 

In the present paper we finish the calculation of the leading nonperturba- 
tive (NP, henceforth) contributions to the n = 2 states. We are thus able to 
present a complete, rigorous and parameter-free QCD evaluation of the full 
n = 1 and n = 2, I = 1, bottomium system. For some of the quantities the 
NP corrections (which are always large) are under control; for some others 
the calculation loses reliability. By and large, nevertheless, a coherent picture 
and good agreement with experiment are obtained. 

NP corrections grow very fast with n so for n > 3 they get so large (for 
bb) that a QCD calculation based on leading effects becomes meaningless as 
was indeed to be expected. However, we present some results for n = 3,4,5 
with a view to future applications to the ti system for which NP corrections 
remain small up to n ~ 5. 

This paper is organized as follows: the perturbative qq hamiltonian is 
reproduced in Sec. 2 for ease of reference. The NP corrections to the interac- 
tion are evaluated in Sec. 3. Sec. 4 contains the ensuing shifts in energies and 
wave functions, which are then applied in Sec. 5 to the complete evaluation 
of n = 1, 2, I = 0, 1, j = 0, 1, 2 and spin s = 0, 1 bound states of bb. The 
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article is finished in Sec. 6 with numerical results and Conclusions. 



2 The perturbative QCD Potential. 



We present here the Hamiltonian for the qq system for ease of reference. We 
write it separating the spin-independent, LS, tensor and hyperfme pieces as 
follows: 



H s .i. 
H (o) 

5 s (/i 2 ) 



^ {0) _ C F (3 a s (/j, 2 ) In rjx 
2n r 
_ C F a s (/j, 2 ) 
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7T 
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7T 
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2 V ; 9 } 18 



D= i(3-*>) + <*-*JiL. 
3 V 2/12 18 

S = Si + ^2 is the total spin, L the orbital angular momentum and 
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n/ is the number of active flavours. The running coupling constant we take 
to two loops, 

2 4vr f _ A lnln/i 2 /A 2 ) 

" sl/i j /3 ln// 2 /A 2 \ /? 2 ln/x 2 /A 2 J ' 

We have lumped the constant piece of the one-loop correction into a s (Eq. @) 
because the ensuing potential is still Coulombic and therefore may still 
be solved exactly. The relativistic, full one loop and leading two loop correc- 
tions to the spin-independent piece are known; see TY for details. We will 
not need them now. The total Hamiltonian is of course 

H P = H S± + V LS + V T + V h{ , (9) 

where the index p emphasizes that only perturbative contributions are taken 
into account. 

A result that we take over from TY is the form of the (spin-independent) 
wave functions pertaining to the Hamiltonian H s i. They are easiest 
obtained with a variational method; one finds that they are given by a formula 
like that for the wave functions of the Coulombic Hamiltonian with the 
replacement of the "Bohr radius", 



by 



b(n,l) 



\n(nid/mC F a s ) + ip(n + I + 1) - 1 
a \ J- H Po^s 



2tt 



, (10) 



*i?(r;a^6). 



A few explicit expressions may be found in Appendix II. 
In particular the wave function at the origin becomes 



*S?(0) 
5 w f (n, 

r (0) 



<°/(0) = {l + 5wf(n ; /)}vl/^(0) 



(0), 



3^o 

47T 



ln( 



mCpa s 



+ ip(n + l + l) -1 



;n) 



As stated, is the solution of the equation 



(C F a s 



An 2 



-m . 



(12) 



When taking into account the full H s - U the energies are shifted to E%\ 



+ ip(n + l + 1) \m 



. , , ■ • , ... • (13) 

47T77/ I mCpds ) 

A last word about the notation: the superindex (0) in say, E^ 
means "of zero order with respect to nonperturbative (NP) effects". 

3 The Nonperturbative Interactions. 

It can be shown (TY and || |4|, ||) that the leading NP interactions, at short 
distances, are those associated with the gluon condensate; and, of these, the 
dominant ones are those where two gluons are attached to the quarks. These 
interactions are equivalent, in the nonrelativistic limit (including first order 
relativistic corrections) to those obtained assuming the quarks to move inside 



a medium of constant, random chromoelectric, £ and chromomagnetic, B 
fields. Because the fields are constant they may be considered to be classical; 
and because they are random we may take them of zero average value 

(£) = (B) = 0. 

The average is taken in the physical vacuum. Quadratic averages are non- 
vanishing and may be related to the gluon condensate. With i, j spatial 
indices and a, b color ones one has (for N c = 3 colors) 

(g*BlBl) = -(g^SZ) = ^^W) . (14) 

The relativistic interaction of a quark (labeled with index 1) with classical 
vector fields may be described by the Dirac Hamiltonian 

H D1 = ia t ■ Vi -07 • A(fi) +/3im, (15) 

N 1 = tt-A^ being gluon fields (in matrix notation). A convenient gauge 

a 

is that in which 

A^-n-S, A 1 = - l -r 1 xB . 

To solve our problem one can apply a Foldy-Wouthuysen transformation p] 
to obtain the Hamiltonian (correct including first order relativistic effects) 

H FW1 = m+ -L(p 1 -( ? 4 1 ) 2 --l J p 1 4 (16) 
-°-S 1 -B--^S 1 -(£xp 1 ), 

Si the spin operator and pi = — zVi. Adding to this the Hamiltonian of the 
antiquark (g — ► — g, fj. — > r* 2 ) and their interactions given in the previous 
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section we find the full hamiltonian, which now includes leading NP effects, 

H = H'">-^^ + V LS + Vr + V M (17) 

-gr-g + 7^~ 2 {S xp)-£~ -{Si ~S 2 )-B. 
2m z m 

H(°\ Vls, Vt, Vhf are given by Eqs. (^) to Some of the peculiarities 
of Eq. (pTT|), in particular the absence of an L • 5 interaction as well as the 
presence of a term involving the differences of the spins, had been noted in 
the similar case of the Zeeman effect in positronium [[jj . In Eq. ( |T7D we have 
omitted a term obtained when expanding the square (pi — gAi) 2 in Eq. flIE|) , 
viz., the piece A\ . It would have produced a term ir(a s G 2 ) r 2 /(A8N c m), to 
be added to Eq. (P). The reason for its omission is that it gives subleading 
corrections to all processes (as compared to the contributions of the other 
terms). 

Before embarking upon detailed calculations, let us elaborate on this mat- 
ter of leading and subleading corrections. Because 

2 



(p) ~ mv ~ itlCfOi 



s 



it follows that the NP terms in Eq. (|17|) are 

£p0xp)'g ~ 5., (18) 

——{Si — S 2 ) ■ B ~ (5.)°. 
m 

This simplifies enormously the calculation at the leading order as seldom 
more than one, and at most two terms need to be considered. A further 
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simplification is that, with the only exception of the hyperfine splitting for 
n = 2, I — 1, only the tree level piece of H p has to be taken into account 
when evaluating leading NP effects. 



4 Energy and Wave Function Shifts. 

4.1 Spin— independent Shifts. 

Although most of the spin-independent shifts of energies and wave functions 
were discussed in TY and [[|, we give here a detailed calculation for ease of 
reference, to correct an error common to TY and to present the results 
for the n = 2 wave functions and to explain in this simple case the way the 
calculation works. 

The effects of the nonzero condensate are evaluated with the help of 
perturbation theory. The perturbation consists of the terms (cf. Eq. (|T 



— gf ■ £, —^t(S X p) ■ £ , — — (Si — 5 2 ) • B. Because, for spin independent 
2m z m 

effects, the first term gives a nonzero result we may neglect the others which 



would contribute corrections of higher order in a S} cf. Eq. (p!8|). Second 
order perturbation theory is required as only quadratic terms in £ will give a 
nonvanishing contribution, as discussed in the previous section, Eq. (|14]) and 
above. The method of evaluation, for this particular case, has been developed 
by Leutwyler, Ref. 0, and is related to Kotani's treatment of the second 
order Stark effect || , up to color and angular momentum complications that 
we now discuss. 
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We denote the solutions of the unperturbed Hamiltonian by 



tf(O) 



(0) \ 
nlM 



= — 

n 9 9 

m or n z 



(o) \ 

nlM I i 

C%a s 



<?m = Y l M {r/r)R^{r). 



(0), 



(19) 



(We have omitted the trivial rest mass energy term). The B$ (r) are identical 

to the standard Coulombic wave functions for the hydrogen atom with the 

2 

replacement of the "Bohr radius" by a = — - ■ Second order perturbation 
theory yields immediately the energy and wave function shifts: 



E = £ (0) +E NP = ¥V + V 



NP 



with 



ij, ab \ 



gnS a t a 



ff (°) - Ei 



(0) 



(0) 
nlM 



(20) 



and 



ij, ab 



ff (°) - M 0) 



if (°) - Pi 0) 



5 ^ ^6 * 



(0) 



(21) 
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Here 



Til 



^(0) \/^(0) 
^ nlM /Y* nlM 



is the projector orthogonal to the nl state. It does not appear in Eq. (^0|) 
because 

/*(°) \r-£\m {0) ) - 



The expressions (p0|) , (|2T|) are first simplified by replacing 

^•••^'--^(« S G 2 ), (22) 

recall Eq. flT|). 

Next we take care of the color algebra. The one-gluon exchange potential 
is given, when acting on arbitrary color states by 

T tti'tk'k ■ (23) 

' a 

If the initial (and final) states are color singlets we may average 



and then we get the potential, and Hamiltonian, 

Cf&s h{0) = _]_ A _ C F a 



s 

! 



r m r 

we have incorporated, as we always do everywhere, the Coulombic piece of 
the one-loop corrections into a s . 

In Eqs. ( PUD and fl2"T|), however, the states V&Sm) are certainly color sin- 
glets: hence the matrices t b (for example) when acting on them will produce 
a color octet state. For a color octet the potential and Hamiltonian are 



O'.s- 



H' {0) = -^ + ^h- (24) 



2iV> ' m 2N c r 
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One then finds 



5 Sabta HW - Ek 0) 



singlet ) = Cp — - 



singlet 



(25) 
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Putting this together with Eq. (^2|) into Eqs. ( ^) and (|2T| ) gives the 
formulas 

7r(a s G 2 ) 



rpNP 



QN, 



(0) 
nlM 



NP 

nlM 



6N r 



P„ 



1 



- Ek 0) ' " 



H'M-EP 

Pnl 



(0) \ 
nlM I i 



(26) 



x £ r - 

% 



i 



(0) 

nlM I i 



(27) 



which takes care of color complications, so we turn to deal with angular 
momentum. Obviously the perturbation is rotationally invariant so the third 
component of angular momentum, M, is not affected by it; but the total 
angular momentum algebra is not entirely trivial. We write 

1 ^ 1 



£ 



(0) 



where A = 0, ±1, and the r\'s are spherical components, 



Using the formulas 



r±i = T-^n ± ^2) , ^0 = ^3 



y^(r» ; ^ = (-1) A Jf Y_ \{r/r) 



and the addition theorem for spherical harmonics we get 



r x Y l M = r £ Cm(/,/',A)F|; +a , 
j'=|/-i|,|/+i| 



Cm(1, I , A) 



2/ + 1 



(Z,M;l,A|f) (Z,0;1,0|Z') 



with 



V 21' + 1 

the standard Clebsch-Gordan coefficients. 



(28) 
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When acting on a function with well-defined angular momentum I we 
have 



1 



(0) 







1 



where 



1 

H = ( r 2 JL* 

mr 2 dr \ dr , 



rr'(O) p(0) 



0. 



/(Z + 1) K« s 



(29) 



(30) 



with k = —Cp for Hj; '' and k = 1/(2 JV C ) for H^. Using this and the 
explicit values of the Clebsch-Gordan coefficients we find that Eqs. (|26|) , 
(|27|) become 

7r(a s G 2 ) 1 



NP 



6iV c 2/ + 1 



x ( i?. 



(o) 

7(i 



R 



NP 



7r(« s G 2 ) 1 „ 1 



i? 



(0) 
nl 



(31) 



x r 



6N r . 2/ + 1 
/ 



P 



nl 



H, 



(0) 



(0) 



p 



111 



+ 



z + 1 



d(0) 



(32) 



We have succeeded in separating the color and angular variables to obtain 
equations involving only the radial variable and radial wave functions. To 
finish the calculations all that is needed is to find the inverses 

1 _™ 1 



H, 



(o) 



E 



R {0) 

(0) • • • 1 l nl > 



H, 



'(0) 



(o) 



R. 



(0) 



This is described in Appendix I. The ensuing expressions for the E^f and R^f 
are collected in Appendix II for a few values of n, I and will be employed later 
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on. The expression we get for E^f agrees with that found by Leutwyler P| 
and also -Rf p , the only wave function calculated in Ref. JJ, agrees with our 
evaluation. 

We have not succeeded in obtaining a closed general formula for B^f (for 
E^f one is given in Ref. ||) but a few general properties may be inferred 
from Eqs. flM]), fl32|) . Because 

= ^(3n 2 -/(/ + l)) ~ ^ 
2 a s 

and each energy denominator yields a factor — — - (see Appendix I) we expect 

a s n 2 

7? NP ^ — r np ^ 77,8 

■"ni ~ ~4 ' ni ~ ~fi ' 

«s a. 

It thus follows that the importance of nonperturbative effects grows very 
rapidly with n. Moreover we expect them to be smaller for energies than 
for wave functions and, generally, to be larger when I = than for I ^ 
(for the same value of n). These properties may be verified explicitly in the 
expressions collected in Appendix II. 

The energies and wave functions correct to leading order in NP effects 
and including one-loop corrections are then 

Enl - E nl + E nl , 

Rni(r) = R^(r) + R7(r), (33) 
*nw = Y l M (r/r)R nl (r), (34) 

the R(°\ m being as given in Eqs. fllQl) , (fTlf) . 
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4.2 Hyperfine splittings 

The hyperfine splittings are caused by the interactions that depend only on 
spin; they are Vm in Eq. (H), and the piece 

--(Sx-S 2 )-B 
m 

in Eq. ([H]). Besides the splitting caused directly by the last term, there is a 
nonperturbative contribution indirectly generated by — gf ■ £ . This contribu- 
tion that we will call "internal", comes about because, when evaluating the 
expectation values 

(*|^hf|*> 

we should use the wave function including the NP corrections discussed in 
the previous Subsection: 

^(^\Vu\^) + 2(^\VmKF). (35) 
The "internal" NP splitting is the last term in Eq. (|35|): 

A™E nls = 2 (¥X\ V M iO • ( 36 ) 
To evaluate this to leading order we use the expression 

V h ^ A -^5(r)S^ 

and thus we get 



A M E n0s = 2s(s + l)^p^R^(0)R^(0). (37) 
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For / 7^ the leading piece of Vm gives zero, because i?^(0) vanishes. We 
have to take into account the radiative correction to Vm and then 



Am E n , 



Is 



2s(s + l) 



AirC F a s ( (3 



3m 2 

,2 D (0)/ 



21 

T 



47T 7T 



o r" 3 



(38) 



It will turn out that, for I * 0, this internal shift will be subleading. This 
fact is very interesting because this is one a the few cases where a rigorous 
QCD analysis yields results qualitatively different from the calculations based 
on phenomenological potentials. This we will discuss in detail elsewhere. 

The contribution to hyperfine splitting of the interation ——(Si — S 2 ) ■ B 

m 

we will call "external" |. It may be calculated as we calculated E^f in the 
previous Subsection. We find, 



A hf E n i s 



[s(s + 1) - 3] 



7r(a s G 2 ) 



p(0) 



H, 



'(0) 



Urn 



d(0) 



(39) 



6N c m 2 

The inverse is obtained with the formulas of Appendix I. 

To the NP contributions we have to add tree level (relativistic) and ra- 
diative ones, that we collectively label perturbative: from TY, 
s(s + l) C 4 F a s {fi 2 )a 3 s (^ 



ALE 



M -^nOs 



X 



2 
1 



3n 3 

2 V mCFa s 



m[l + 5 w f(n,0)Y 



n i 

Y- 



21 



In 



r? 



C pO-t 



-Y-- 



AP J? 

t-^hi 1 -'nls 



8(8+1) 



2 67rn 3 /(/ + 1)(2Z + 1) 



1 

?2—l 

2n 
, 2 



n 



7e 



2n 



(40) 



21 

4 



/ ^ 0. 



§In the case of hyperfine splittings the internal contribution is chromoelectric and the 
external one chromomagnetic, but this is not true in other splittings. 
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The constants are as in Eq.(|]). The full splitting is thus 

Am E nis = A M E n i s + Ajj E rds + E nis , (41) 
with the various pieces given in Eqs. ( PB[ ) to 
4.3 Fine splittings 

Also here we have "internal" and "external" contributions. The internal ones 
are , as before, induced by the NP modification of the wave function. The 
calculation is somewhat complicated because now two operators, the LS and 
Tensor ones (Eqs. (|) and (0) ) contribute. We fincQ 

A j E n ij = 2<W(rM) {(^U + (V^nij} , (42) 

where 

/d(0)| -31 dNP\ 

Mn,0 = U ; (43) 

R^f is given in Eq. ([32]) and Vls , ^ / t°' ) are ^ ne leading (tree level) pieces of 
Vls, Vt- Using the explicit expressions for these we have, 

(vSu - ro + i)-'0 + i)-2l 16n3; S°g + ir ; (44) 
«^ - <5 s ^8^§TT) mi < 46) 

with 

f 7 = I- 1 

1 2Z-1 i J 1 1 

kSWii = +1 ,3 = 1 (46) 
I -27T3 ,3 = 1 + 1- 



^We consider that the states correspond to total spin s = 1. For s = 0, A™E£fj° = 0. 
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The leading "external" fine structure shift, A^ x E n ij, is caused by the 
crossed combination of the perturbations 



-gr-S , 



9 



2m 2 



(Sxp)-S. 



In this case the external shift is also chromoelectric; the chromomagnetic per- 
turbation — — (Si — S2) -B does not contribute to the fine structure. The color 
m 

algebra is now like the one for the spin-independent shift, Subsection 4.1. 
Thus, 

n(a s G 2 ) 



-2- 



X 



2m 2 



(0) 
nlj 



(Sxp)i 



1 



H'W-Es 



(o) 



(0) 
nlj 



• (47) 



The angular momentum algebra, on the other hand, is somewhat compli- 
cated. It is developed in detail in Appendix III for n — 2, I — 1. One 

gets 



AfE 21j 



n(a s G 2 } fj (j + l)-4 / g) 



QN r m 2 



d(0) 
K 21 



(0) 
21 



d_ ( 1 1 



d(0) 
^21 



I„(0) = 1/(1) = -vQ) = ± . 

The calculation is finished using the inverses of Appendix I. The result is 

1780 + 1)- 4]- 2784 z/(j) n(a s G 2 ) 



AfE 2 ij 



K{j) 



9945 
n(a s G 2 ) 
m 3 (CFtt s ) 2 



m 3 (CFa s ) 2 



(48) 
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with 



8976 1 
^(«) = -9945.^(1) = "^(2) = 2^(0). 



(49) 



The perturbative fine splitting is (for s — 1; the splitting should be considered 
to vanish for s — 0) 

3C£ as ( M 2) 5 3 (M 2) 



A? 3 



"0 



+ 



16,«( ; + l)(2 ; + l) ro ^ + 1) -'" + 1) - 2||1+ *- (n ' 0)1 
x jl + (y ~ 2 ) ( lnn - 1 - + ^ + !) + ^(2/ + 3) + ^(2Z) 
n _/_l /2 \ 125-lOn/ ft /i „ J 

C 4 F a s (^)a 3 s (^) 1 2 
8n>i(i + l)(2/ + l) m ( 2 Sl2)y [1 + 5wf (n ' ° )] 

x f! (y -3j (Inn- V»(n + Z + l)+V(2i + 3)+V(20 



n " / - 1/2 ^ ' ^ln-^ + SlnC^ 



n 



7T 



J 2 ITlCfOIs 

The constants as in Eqs. (H), (^) and (p~T|) . 

The full, relativistic plus radiative plus NP fine splitting is then 



Af E n ij — Af i? n y + Af n i£ n jj + Af X £7 n y 



(51) 



the various terms given in Eqs. (f42|), ( f48| ) and (|50|). 



4.4 Decays into e + e 



For a state with Z = the decay rate into e + e is given by 

2 



T{T{nS) -> e+e") = — 



M(T(n5')) 



[mC F a s (/i 2 )] 3 



x (1 + 5 r ) [1 + 6 wi (n, 0) + Pnp(^)] 2 • (52) 
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Here 5 r is a "hard" radiative correction ||, 

AC F a s 



71 



(53) 



5 w f(n, 0) is given in Eq. ( |TTD and Pnp(^) is the ratio of NP to unperturbed 
wave functions at the origin: 

g(Q) 

It is to be calculated with the expressions of Appendix II. 



PnpW = -fcf • (54) 



5 Properties of Bottomium in States with n = 

1, 2. 

We will use spectroscopic notation: states will be labeled n 2s+1 lj, / = 0, 1, 2 . . 
or S, P, D, .... The somewhat whimsical notation of the Particle Data Ta- 
bles (PDT) [|TI]] will also be indicated. For n = 1, 2, 3 mixing does not 
occur. 

5.1 States with n = 1. 

From TY we have 

111(1%) = M(T) = 2m <j 1 - 



1 Cj,5 a V) C 2 F (3 a 2 s (^)a s (^) 



8 87T 

x r^sr +1 - 7E J) + 7s^F m - (55) 

1872 

eio = ~ 1.468 . 

10 1275 
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The order af is partially known|J it adds to the right-hand-side of Eq. (|55|) 
a term 



2m 



3C F 
16 



1 + ai + 



7eA)\ « 



— a s a s + 

7T / 



2~2 



5C^ _ 4 
-a„ — 



In 



- 1 - 7E ) a ^ + it a ^ 



.(56) 



128 ~" s 16tt 2 

We will use both Eq. (|55|) alone and Eqs. (|55|) plus (p6|). 

The hyperfine splitting is obtained from Eq. (|4~T|), (Eq. (|37D) evaluated 
with the expressions for the i?'s of Appendix II, and the inverse in Eq. (|39|) 
with those of Appendix I. The result is 



M(1 3 S!) -Mil 1 So) 



M(T) - M( Vb ) 
C F a s (fi 2 )a 3 s (fi 2 ) 



■ m 



x 1 + 



A) 



In 



\i \ _ 21 / 1 

ITiCfOLs 4 V CfO-s 



C F a s {fi 2 )al{fi 2 ) m 



270459 1161 



108800 8704 



-lj+B 

n(a s G 2 ) 



71 



(57) 



In the NP contribution the first term is the internal, the second the external 
which is, as is generally the case, substantially smaller than the first. The 
term in square brackets is, after multiplying by 7r, 7.81 + 0.42 = 8.23, slightly 
smaller than the value given by Leutwyler || which was also used in TY and 
equal to 10.2. The difference in the value of the hyperfine splitting, however, 
is fairly small. The corrected value, following from Eq. (|5"7j) , will be given 
below. 



"it includes leading relativistic corrections 0(aj.) 7 

one-loop radiative ones 0((a^/ir) In// 2 ) and C(a^/7r), and leading logarithm two-loop 
corrections 0((a 4 s /tt 2 ) In 2 ^ 2 ). The error of Eq. <E% should be at the 10 to 20 % level. 
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For the e + e decay Eq. (|52| ) gives us 



r(l 3 Si ^e + e" 



T(T -> eV 



M(T) 



mCj?a s (// 



7T 



X 



1 + 3/? In 



2/x 



a s 270459 n(a s G 2 ) 



mC F a s ' 2 ,J 7 4tt 217600 m 4 5f 
and we have inserted the explicit values for S r , <5 w f, pnp- 



(5? 



5.2 States with n = 2. Spin— independent shifts. Decay 
into e + e . 

We will denote by M(2 3 P) the average of the masses of the states** 2 3 Pj, j = 
0, 1, 2: 

M(2 3 P) = - |5M(2 3 P 2 ) + 3M(2 3 Pl) + M(2 3 P )} = 9900 ± lMeV . 

(59) 

From the analysis of TY and Ref. || we have, 



M{2 3 S 1 ) - M{l 3 S l 



M(2 3 P) -M(l 3 Si 



M(T(2£)) - M(T(15)) 
r 3C 2 5 s V) 



2m 



[ 32 

32 



3 In— ^ + --3 TE -m2 



7T 



(2 6 6 2 o-e 10 )7r(a 5 G 2 ) 2102 
+™ — A — — , e 20 = ~ 1.585 ;(60) 



J F 

\ 32 

C F (3 a s a s 
32 



1326 



2m 



3 In „ - - + — - 37e - In 2 



CF?Tia s 6 



7T 



'* Denoted by Xbj0-P) by the PDT people, Ref. M 
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(2 6 e 2 i - e 10 )7r(a s G 2 } 9929 , , 

+m - „ A /J — - , €21 = 7^ ^ 0.9984 (.61) 



9945 



It is interesting to consider on its own the " Lamb shift" , difference between 
Eqs. ( |6"0"D and (|6T|), as here only the states with n = 2 are involved: 

2 6 {e 20 -e 21 )ir{a s G 2 } 



M(2 3 S0 - M(2 3 P) = 2m C P^ a ' S ' 



+ m ■ 



96tt 

As for the decay T(2S) — > e + e~, Eq. ( |52"D gives 

2 

K F a s (/i 2 )] 3 



Clm 4 5 4 



(62) 



r(2 3 5 



x 1 



l^-e'e ) 
AC F a, 



7T 



M(T(2S)) 
2/x 



1 + 3/3 In 



mCVa s 2 



1 \ a s 302859 7r(a s G 2 >TL 
^- 7E j47 + ^84-^^ (63) 



5.3 States with n = 2. Fine splittings. 

From Eq. ( |5"ID and after some work we get the fine structure splittings^ 



M(2' i P j ) - M(2' A P) = m^ a s (|i 2 )^(/i 2 ) 



x 



X 



2(JL 



l + 3/5 In . , 
V mGpOt s 6 

j(j + l)-4 



7e 



+ 



256 
384 



1 + 
1 + 



£-2 
2 



47T 

In 



1" 

2// 



111699 %{a s G 2 )' 



221 



'A 



2\i 



, n . fi 125 -10n/ 
7e +21n— + - 

m 



36 



7T 



In ?~„ +l-7 E )+31n^ + P 





1] 


-J 





+ m 



ir(j)7r(a s G 2 ) 



m 4 (CV5 s ) 6 

The first term containing (a s G 2 ) is the "internal" NP shift (corresponding to 
Eq. fl42])); the last term is the "external" piece, Eq. (|48|). The experimental 



(64) 



^Because <5np, <5 w f are large we have included them in a factor [1 + S w f] 2 (1 + 2<$np) in 
Eq. (p4]). This form or the equivalent one of a factor [1 + 5 W { + <5np] 2 are the ones that 
give more stable numerical results. 
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shifts are 

Af (2 3 P 2 ) - M(2 3 Px) = 21 ± 1 MeV 
M(2 3 Pi) - M(2 3 P ) = 32 ±2 MeV. 

5.4 Hyperfine splittings for states with n = 2, / = 1. 

The hyperfine splitting M(2 3 P) — M(2 *Pi) has not been measured experi- 
mentally for bottomium. For charmonium, 

M cg (2 3 P) - M cc -(2 ^ = -0.9 ± 0.2 MeV . (65) 

The theoretical calculation has been displayed in Subsection 4.2. After sub- 
stituing the explicit expressions for the various pieces we get 

'0a 21 \ C^alal 



M(2 3 P) -M(2 x Pi) 



m 



2 4 J 288tt 
6l7r(a s G 2 ) 
117 m 4 al 



+ m 11,^2 • ( 66 ) 



0o 21 

This effect is remarkable. The coefficient — — is negative; hence the per- 

turbative and all internal NP contributions (which are, however, subleading) 
will be negative. On the other hand, the external NP correction is positive. 
For the (relatively) light quarks cc, the perturbative piece dominates; but for 
bb, because it decreases like af, and the NP one grows like a~ 2 , the situation 
is reversed and we will get 

M 65 (2 3 P) -M 65 (2 1 P 1 ) >0. 

This is of importance for calculations based on phenomenological potentials 
(see e. g. Refs P, PD, a matter that will be discussed in a separate publica- 
tion. 

24 



6 Numerical Results. 



The numerical results which correspond to the formulas given in the previ- 
ous sections are presented in Table. I. Before discussing them a few words 
have to be said about the calculational procedure. The quantities pertaining 
exclusively to bb in states with n — 1 have been taken from TY with the only 
modification of the hyperfine T — i] b mass difference where we have incor- 
porated the (minute) modification following our corrected evaluation of the 
NP contribution. The criterion adopted in TY to choose the renormalization 
point /i, was to require that radiative and NP contributions be equal in ab- 
solute value. Most results were in fact little dependent on the actual value 
of \i chosen. The reason is that, for n — 1 the quark mass (as a function of 
M(T) taken as input) begins at order a° s and the first corrections are 0(a 2 s ). 
For the decay T — > e + e~, the leading contribution is order a 3 s ; finally the 
"Balmer" mass differences M(TnS) — M(T1S) start at order a 2 . By con- 
trast the Lamb shift M(2 3 S 1 ) -M(2 3 P) starts at 0(a 3 s ), the fine splittings 
among 2 3 Pj states begin at order at. (as does the n = 1 hyperfine splitting) 
and, finally, the hyperfine splitting M(2 3 P) - M(2 1 P 1 ) is an effect of C(a 5 s ). 
This means that for all these quantities the choice of \x is essential as small 
variations in \i get amplified. Because of this we have chosen to fit the value 
of \i. We have considered three possibilities: fit the two fine splittings, and 
then the Lamb shift and Balmer splitting M(2 3 5i) — M(l 3 S'i) come out as 
predictions; include the Lamb effect in the fit; or fit all four processes. We 
present results in all three cases; we consider the last possibility to give the 
optimum calculation. A remarkable fact that lends credence to our results 
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is that the values of fi obtained with the three methods, as well as with the 
criterion of TY (for the Lamb shift and Balmer splitting that was considered 
also there) are extremely close one to the other. 

The values of A, (a s G 2 ) were not fitted. We chose, as already mentioned, 

A(n/ = 3, 2 loops) = 250 j^MeV, 

(a s G 2 ) = 0.042 ± 0.020 GeV 4 . (67) 

Because we take M(T) as input, we deduce (and fh^fh 2 )). For the pole 
mass, Eq. (|67|) implies, according to the analysis in TY, 

m b = 4906 tll(A) +t((a s G 2 )) MeV , (68) 

the first variation in Eq. fl68|) tied to the variation of A in Eq. (|67D, the second 
tied to that of the gluon condensate also in Eq. (|B"7|). 

The agreement between theory and experimental data is remarkable, as 
is remarkable the stability of the predictions of the (as yet unmeasured) 
hyperfine splittings. The deviations are of the expected order of the higher 
corrections, 0(a s ) ~ 30%. As drawbacks, however, let us mention the fact 
that some of the NP corrections, notably the ratio <5np, do actually exceed 
unityQ. This makes the results of the fine splittings less impressive than 
what they look at first sight. Nevertheless, the choice of [i as well as the way 
to write our equations certainly allow a control of the results. 

The process T(2 3 S*i) — > e + e~ merits a special discussion. If we take the 

central value /i = 976 MeV (Table I, column(c)) and consider the leading 
"A list of some radiative and NP contributions is given in Table II. 
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expression of the width, i.e., we neglect radiative and NP corrections, we get 



r (o) = I 
4 



Q b a 



2 



C d F m d ai = 0.77 keV . (69) 



M(2S) 

This is the value reported in Table I, and it compares favorably with exper- 
iment. Unfortunately the corrections involve the factors 

(1+S r ), (1 + M2,0)) 2 , (1 + Pnp(2)) 2 

(see Eq. fl63f) for the expressions for the 5, p) and one has 

8 r = —0.61 , <5 w f = —0.53 , pnp — 3.6 . 

The prediction then becomes meaningless since the corrections are much 
larger than the nominally leading term, Eq. (|69"D ; although here, as it happens 
in the cc case (see TY) this leading term yields a reasonable evaluation, 
considered as an order of magnitude estimate. 

Taken all together, our results here as well as those of TY, constitute a 
coherent description of the lowest lying states of heavy quark systems, using 
only rigorously derived QCD properties and without need to have recourse 
to phenomenological potentials or adjustable parameters. 
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Appendix I. 

We evaluate the inverses 

1 



Here 



2^ „, 0) 1 C 2 a 2 



na ' ma 2 n 2 An 2 mCpa a ' 
and 

RK = 1 1 d ( 2 d\ | 1(1 + 1) | K g s 

mr 2 dr \ dr J mr 2 r 

For v integer p u turns out to be a polynomial: 

u+l 

and 

Cf mn 2 a 2 



Kn + (u + 2)C F 4 ' 



9-i = [j(j + l)~ 1(1 + 1)} Cj , j=l, l+l, v+1; 

Cj = , j < I . 

'(0) 

When Hi = H l ' those equations give a unique well-defined p u . For Hf = 
H^ one should replace n by n + e. Then p u contains a singular coefficient, 
proportional to 1/e. However, when evaluating 



H^ - nl P 

with P n i the projector orthogonal to the solution of 

(hp - £«»>) <» = o , 

the singular term drops out and the limit e — > may be taken. 
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Appendix II. 

Here we list some nonperturbative energy shifts and wave functions (spin- 
independent). We write 



tti„i — — ; — „ ~ ; : — m 



e n ; n b 7i(a s G 2 ) 



J nl 



Then, 



624 1051 

e 10 — "TTTr 7 e 20 



425 663 

9 929 769 456 

621 ~ 9 945 C3 ° ~ 463 239 

11562 272 101509 

C 31 — o ^ r QO 

e 50 



8 492 715 u 60 060 

443 288 368 



260175 675 

2r 

For the wave functions, and with p = — . 

na' 

pNP _ 7r(q s G 2 ) 2 /2 



m 4 (C F 5 S ) 6 a 3 / 2 



2 968 104 2 52 o 1 4 

p 2 p 6 p 4 

425 425 1 1275 H 22W 



NP n(a s G 2 ) 1 /2 

3 828 736 1914 368 134 528 2 67264 „ 736 4 16 = 

p p H p H p H p 

1989 1989 H 1989 F 5 967 ^ 663 H 153 F 

oNP = 7F (^ G ' 2 ) 1 „„- P /2 

r 3 299 840 149 888 2 5 248 3 32 4 
X I 1 989 5 967 P ~ 1 989 P ~ 153 P 
pNP _ K(a s G 2 ) 1 /2 



m 4 (C F 5 s ) 6 ^a 3 / 2 

f 189 965 808 189 965 808 24 735 864 2 3 462 552 3 1302 4 
X l 5719 5719 P+ 5 719 P + 5 719 9 ~^3~ P 
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3 042 5 9 6 
P To P 



1 505 r 43 

dNP k{(X s G 2 ) 1 p /2 

31 ~ m 4 (C F 5 s ) 6 v^ a 3/2 



f 1325 287 104 331321776 2 124 833 216 , 49 872 4 

x <^ p p l p 3 H p 4 

I 62 909 H 62 909 ^ 314 545 H 1 505 ^ 

3 672 - 9 f 



dNP Tr^G 2 ) 1 -p/2 

40 m 4 (C F 5 S ) 6 a 3 / 2 

f 5 609 365 504 2 804 682 752 57 706496 2 20160 512 

x <^ p H p z p 

I 45 045 15 015 H 1001 F 15 015 H 

93 551104 , 59 392 - 256 fi 32 7 

p H p H p H p 7 

135135 r 3 861 P 429 r 35 r 

nNP _ n(a s G 2 ) 1 /2 

-n-50 — : ; ~ ~ — i= _ - e 



3 



m 4 (CV5 s ) 6 v^a 3 / 2 

37087558150 000 74175116 300 000 35 702 282 000 000 
+ oT^TTT^ oT^TTT^ P 



31221081 31221081 r 31221081 

13 695 312 550 000 „ 561983 427 500 „ 138 527387 500 



P- no^ooxo P + 



93 663 243 ' 93 663 243 ' 93 663 243 

4 827500 „ 1250 7 625 8 



261873 P 9 699 P 6 588 P 
For ease of reference we also give the first R^'s 



(V 



4?(r) = * - -e-/ 2a 
21 V ; 2 v / 6a 3 / 2 a 



L 20 



v/2a 3 / 2 



a 



For the -R^'s, replace a by 6(n, Z) given in Eq. (|10]) . 
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Appendix III. 

We evaluate the matrix element (21 stands for nl) 



(o) 

21j 



(SxP) l - 



i #'(o) - Ei 0) 



(0) 
21j 



It is convenient to use a Cartesian basis for the spin-angular momentum 
piece of \I>2ij so that 



(0), 



ik 



Here r — f/r, the \k are column spin 1 wave functions and the coefficients 
are 

d 0) (o) = -L^ , d a) (i) = ^^ , 

V47T V0 7T 



df(2) 



'An I 3 



The last expression valid for a ^ b. The indices 0, a, ab, collectively de- 
noted by a in (III. 1) give the (Cartesian) third component of total angular 
momentum. The spin-angular momentum wave functions 

(<*), 



f(j) = E41j) hxu 



ik 



form an orthonormal set: 



(III. 1) 



dsu (a) 0'K w 0") = s^s^. 



We have 



A* 



a \ 
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ik i'k' 
abc 



X 



r a hiik{j)Xk 



1 



d(0) 
K 21 



If we write identically 



[fafi - \ 5 a?j + ^ Sai , 



then the first term corresponds to angular momentum 2, and the second to 
angular momentum zero. Therefore, when acting on the first we may replace 
if by H 2 < ~°\ and when acting on the second H'^ by H ^°\ Hence, 



M = £ (R 



(0) 
21 



iki'k' 
abc 



J dVt & k , r v x\> tabc Sb P c 



X ( r a fi ^ 5 a i J ^'(0) _ ^-,(0) ^ikXkf 



d(0) 
K 21 



+ 5 E <* 

aftc 



iik XkT 



d(0) 
K 21 



and, after straightforward substitutions and arrangements, 



M 



E (a 



(0) 
21 



iki'k 1 
1 



iik Xkr 



d(0) 
-^21 



iki'k 1 cs \ 



x 



9_ / 1 1 

dr Uo' (0) -4 0) /f 2 ' (0) -i?f, 



d(0) 
^21 
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'(0) 

The only noteworthy aspects of the derivation are first, that, because H l 
only acts on the radial variable, and the fj only depend on the angular ones, 

1 . „ 1 



and, second, that for any /(r), 



(0) P (0) ' 



Pkf{r) = -if k 



dfjr) 
dr 



The calculation is readily finished. Because 



J2 tikU) ?iXk 
ik 

corresponds to total angular momentum j, 

? ? rt/v j(j + l) -1(1 + 1) -s(s + l) x , ( 



ik 



ik 



with I = s = 1. Defining also 

4tt 



9 



(^k(MkiU) - Cu(Mkk(j)) , 



\v(0) = 1/(1) = -1/(2) = 1 , 



we finally get 



X = 4 ~ i( i + 1) /iff 



r tf 2 ' (0) - £f 



d(0) 
^21 



(0) 
21 



9 



<9r Vi/ ' (0) -4 0) 



p(0) 
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Quantity 


(a) 


(b) 


(c) 


Experiment 


[l (MeV) 


qqn + 213 " 43 

yyu -198 +90 


Qfi o +231 -53 
yOo _224 +104 


Q7f y +238 -54 
y '° -228 +107 




a s 2 ) 


0.36 ±0.03 t°o°ol 


0.37 ±0.02 t° ° l 


0.36 ±£^±0.02 




a s {\x 2 ) 


n c/l +0.07 +0.02 
u - d ^ -0.06 -0.04 


net; +0.06 +0.03 
u -°° -0.03 -0.05 


n ^ +0-05 +0.03 
u -°° -0.03 -0.05 




2 3 P 2 - 2 3 Pi 


22.2 -° A 2 ± 


-i o Q +2.3 -2.6 
lO.y _6.1 +2.5 


on n + 1 - 2 _2 - 5 

zu - u -6.7 +2.6 


21 ± 1 MeV 


2 3 Pi _ 2 3p o 


30.0 ±0.6 ±8.2 


zd -° -9.0 +3.1 


97 9 +5.2 —3.6 
-9.8 +3.2 


32 ± 2 MeV 


2 3 Si -2 3 P 


1 no -49 +92 

iyo +82 -so 


183^40 tH 


ioc -39 +32 
loD +40 -42 


123 ± 1 MeV 


2 3 Si - l 3 ^! 


/1 87 ~ 148 + 68 
+222 -69 


Ana -105 +14 
430 +74 _28 


4 cc -97 +17 
+68 -32 


563 ± 0.4 MeV 


2 3 P - 2 ^ 


i.7 +8:? ± 0.7 


1.6 +g;t ± 0.6 


1 7 -0.6 +0.5 x, 
1./ +o.4 -0.6 MeV 




2 3 S 1 -*e+e- 






~ 0.77 


0.56 ±0.10 MeV 


m b (m b 2 ) 


4397 t£ ;1 MeV (d) 


4250 ± 100 (f ) 


l 3 ^i-l% 


33 l* 3 t\ MeV (e) 




1 3 5 1 -> e+e- 


1.01 ±0.02 toilkeV 


1.34 ±0.04 


1 3 S 1 ^2 7 


0.17 keV (d) 





Table 1: Compilation of results. 

Theoretical predictions, and experimental values for bb states with n — 2, 1 
and I = 1, 0, s = 1, 0, j = 0, 1, 2. 

(a) The parameter jj, obtained by fitting 2 3 Pj. 

(b) Fit including also 2 3 Si - 2JP. 

(c) Fit with the former and 2 3 Si — 1 3 S\. 
(x 2 /degrees of freedom = (0.33 +° f 2 +° g) /S). 

(d) Result from TY. 

(e) Result with analysis from TY with corrected NP contribution (N. B.: 
old result, 35 MeV). 36 

(f) Values obtained from e + e~ — > hadrons via QCD sum rules, see Refs. 4. 



Quantity 


tree (a) 


tree + rad. (b) 


NP ext. (c) 




5np 


Total 


2 3p 2 _ 2 3 p x 


11.5 


3.4 


1.9 


-0.27 


2.2 


20.0 


2 3 Pi _ 2 3p Q 


14.4 


4.9 


0.95 


-0.27 


2.2 


27.0 



Table 2: Sample set of contributions. 

with n = 976 MeV ; A(n f = 3, 2 loops) = 250 MeV ; (a s G 2 ) = 0.042 GeV 4 . 

(a) With tree level potential (including relativistic corrections). 

(b) One loop radiative corrections. 

(c) External NP corrections. 
All dimensional numbers in MeV. 
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